NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL MEMORANDUM 

No. 1162 

TUNNEL CORRECTION FOR COMPRESSIBLE SUBSONIC FLOW 

By A. V. Baranoff 


Translation 

“Zur Frage der Kanalkorrektur bei kompressibler Unterschallstromimg” 
Deutsche Luftfahrtforschung, Forschungsbericht Nr. 1272 



Washington 
July 1947 



KAnONAL ADVTSORy COMMITTEE FOR AERORAOTICS 


TECHNICAL MEMfMHDUM NO, ll62 

TUMEL CCREECTION FCE COMPRESSIBLE SUBSONIC FLOW* 
By A. V, Barcaioff 

SUMMARY 


This report presents a treatment of the effects of the tiainel 
walls on the flow velocity and direction in a canpressihle medium 
at suhsonic speed hy an approximate method. Solutions with numerical 
calcxilations are given for the rotationally symmetric and two- 
dimensional proIxLems of the flow past hodies, as "vrell as for the 
downwash effect in the tunnel with circular cross section. 


1. SYMBOLSi 

h wing span of the model wing 
r circulation 

h half of the tunnel height, two-dimensional case 
^ profile volume of the model, two-dimensional case 
p tfcich number sq,uared in the undistur-bed flow 
q. variable of integration 

R tunnel radius 

o variable of integration 

T volume of the model, case of rotational symmetry 

*"Zur Frage der KanaUcorrektur bei kcmpressibler Unterschall- 
stromung, FB 1272, Zentrale fur wissenschaftlic.hes Berichtswesen 
uber Luftfehrtforschung (ZW3) Berlin-Adlersiiof, July 19^. 

^This list only contains . symbols appear’ing In the final results 
(©q.uations (17), (25), (31), (32), (4l), inad (42)), Symbols used 
in intermediate calculations are explained at the point of their 
introduction. 


2 


MCA IM No. 1162 


U flow velocity 

u increased velocity at the tunnel wall 

u* additional axial velocity (due to the constriction of the stream) 

w* additional upwash velocity (due to the constriction of the stream) 
^ , p, or § , T] are the coordinates for the case of 
rotational symmetry or two dimensions; rendered di.mensionless 
"by division "by E or h. 

2. GENERAL STATI3‘IE3!W’ OF TBE PROBLEM 


The effect of the tunnel walls on the flow around a "body 
acq,uires increased significance at high velocities as much throvigh 
compressihility as through the often unfavorable ratio of model 
dimensions to tunnel diameter. In this, the question concerns 
the effects on the flow speed and direction, the first case of 
which is, possibly, that of a model, symmetrically suspended, in 
a flow where there is zero lift, while the second case Is that 
of a circulatory flow past a thin profile. The differential 
equation for compressible subsonic flow should be taken as a basis, 
here, in the approximation form named after Prandtl. If 
represents the velocity potential, in cylindrical coordinates this 
equation, then, reads: 


5^ 

Sr^ 



1 ^ 2 $ 
r-2 


+ (l - n) 


5x2' 


= 0 


( 1 ) 


This holds for a so-called, near* parallel flow, that is a 
tmifona principal flow in the direction of the x-axis on which 
is superimposed a flow of ordinarily small veloci1;y'. 

Now let be the potential of the flow in the medium, 
unconfined, and ^ the potential of the additional flow appearing 
because of the effects of the tunnel walls . I’**’ certainis'- satisfies 

the differential equation (l) in the entire range of the interior 
of the tunnel as a good approximation. The same cannot be said 
of $ because in the vicinity of the body the deviations from the 
principal flow can be of the same order of magnitude as the 
principal flow, itself. The quantity, §, will, however, certainly 
do at a distance from the body, that is, in the neighborhood of 
the tunnel wall, possibly, ;just as well as * of equation (l) . 
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Now since the connection "between $ and consists of the 
fact that 


-|- ($+$*)= 0 ( 2 ) 

at the tunnel wall, the desired is touched only slightly 

hy the uncertainty in the potential ^ in the vicinity of the 
"body as far as it succeeds, that is, in giving solutions of 
equation (l) of such a kind which describe the action of the body, 
which the flow moves past at a great distance from it with 
sufficient accuracy. First of all, in the following the rotationally 
symmetrical and the two-dimensional problem’ for the flow past a 
model will be treated, for which ascertaining a correction factor 
for the flow velocity or its Mach nxaaber is the object of this 
investigation. In the conclusion, the problem of downwash correction 
factor is handled in connection with that. In a formal sense the 
method in all three cases depends on the same ai’tifice (compare 
reference l) , namely, in that the condition, at the edge (reference 2 ) 
is satisfied, first of all, within a finite longitudinal section 21 
of the tunnel cylinder, and the limit Z— is taken only then. 

The solutions all appear, therefore, in the form of Fourier integrals. 
It should be mentioned that in the two-dimensional case the method 
of reflection of the singularities (reference 2) leads to a solution 
that is more convenient for the purpose of numerical calcxlLation. 


3 . EFFECT OF THE TUNNEL WALL ON THE FLOW PAST BODIES 
(CASE OF EOTATIONAL SYMMETRY) 


It is logical to describe the disturbances that a body past 
which there is flow, causes at some distance from itself by a 
superposition of sources and sinks in which the sounce and sink 
potentials satisfying equation (l) are readily expressible# The 
discussion is limited, at this point, to the case where the body 
is s mall enough in comparison to the tunnel radius so that its 
action can be replaced accm-ately enough by that of a sirigle 
dipole. The potential of such a dipole with the x-axis as its 
axis of symmetry, figure 1, reads: 


,®r 

l/'x2 + (1 - u)r2 



k 
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Eegarding the meaning of the dipole moment aij,, arguments 
will not he presented till section 5* 

For the additional potential $*, which gives the action 
of the tunnel walls on the flow, the followi3ig estimate is made: 

~ P(r) X(x) 

with which the following equation derives from (l) 

P” + ip» + (1 - n) P ~=,0 (h) 


First of all, to satisfy che boundary condition (2) only 
for |x)<2, it is necessary to set 


X" + X = 0 


( 5 ) 


in which 



k = 1,2,3 


It is readily seen, that because of (3) and (2) x appears to 
an odd power in so that only 



( 6 ) 


enters in as a solution of (5)* On account of (5) equation (4) 
transforms to 


P" 



(7) 
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The solution of this so-c&llei' Modified- BessejL^ s differential 
eqLuation is 

: ( 8 ) 

where Iq is the modified Bessel function of the first type and 
zero order. The correspondinf» function of the second type does 
not enter into the question hecauBe of the requireiaent of 
regularity for The general solution develops from (6) and 

(8) hy summation over all integral values of' fc. With the use 
of the dimensionless quantities 






(9) 


it reads 



(10) 


The definition of the coefficients follow from the 
"boundary condition (2) . To begin with, for p = 1 


~ krt ^ kJt \ ^ kng 3/1-1^ 


Bn 




+ 1 - u 


Expanding the right-hand side in a Fourier series in by 

comparison of coefficients, after seme intermediate calculation, 
the following is obtained 

X kjta 
cos -r — da 
K 



"k 


vi * ii DU ir~P 

^ Va^^ + l- u 


(12) 


2«2r2 


'o' 
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Now sutstituting (12) in (lO) and talcing the limit X->oo the 
following is obtained 


^ 

0 


sin (q,^) Ifi ( \/l - |i q.P)<iQ cos(q.cx) iei 


Iq’ ( VI - u a) 


^ 0 \j aP' + 1 - [I ^ 


( 13 ) 


The inner integral in this can be put in a form where it is 
expressible by a modified Bessel* s function of the second type 
and first order. This is, namely, 


f cos(qa) da qK3_( VT - U cl) 

° Va^ + l- u^ \ll ~ \i 


where K^_ is the fvinction mentioned. Equation (13) reduces, 
by means of this, to 

/o® ^ 

m^ / sin(q§ Iq ( \ll ~ U <ip) T % ( \Jl ~ u q) dq 
$ * = — ~ / (14) 


? p 

2n"R‘=^ 


I t ( \Jl - [I q) 


A new variable of integration can be written for \/l - n q 
here (represented again by q in the followii^L;) and the expression 

%■ r / »s\- , V 

U* = ■ ' n / COB I “7==r]I (qp) 

2rt%3^1 . \]/r- J ° I (q) 


dq 


(15) 


is obtained for the axial additional: velocity. In 'vri.nd tunnels, 
there is the possibility of finding the velocity at the tunnel wall 
by measurement of the static pressui’e at the tvamel wall. The 
increased velocity there is computed from the two potentials (3) 
and (l4) for ^ = 0 as 


u = 


“r 


1/1 - u ' 


n 

— + 



q%i(q) 
In (q) 


( 16 ) 
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Ulimlnating the dipole mcsnent in ( 15 ) with the aid of (I6), 
then at the position of the tody = P = O) 

Uq* = 0*45^»-u ( 17 ) 


is obtained for the con*ectlonHractor velocity where it can he 
ascertedned by measurement. (See section 7 for the numerical 
calcxilation of the factor, ) The relationship (l 7 ) is independent 
of Mach number. 


k: EFFECT OF THE TUNNEL WALL ON THE FLOW F.AGT BODIES 

(Two -Dimensional Case, see Fig. 2 ) 

In the two-dimensional case the differential eq.uation for the 
velocity potential reads 


d 2 $ ^ 2 $ 

— ^ (1 -P) ^ = 0 


(18) 


A solution of this eq,uation, which is associated with the 
dipole, reads 

ma X 

$ = - 

2n + il - M-)y'^ 


First of all, the moment m^ is simply regarded as given; later on, 
its relation to the size of the body and to the Mach number will 
be discussed further. To fvilfiiJ. the boundary conditions at the 
upper and lower tunnel wall (y = th) req.uires the introduction 
of an additional potential which likewise should satisfy 

equation (18) . 

Through the statement 


-y(y) X(x) 
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it is easy to get a general solution. Its form consistent with (19) 
and the houndary condition reads 


q)* cjj cosh^\/] 


k7ty\ 

I 


sin 


krtx 


( 20 ) 


To satisfy the houndary conditions, exactly the same pi'ocedure 
is to he’ ohserved as in the preceding section. After taking the 
limit as Z->w and hy applying the dimensionless formulas 



( 21 ) 


the expression. 

* - ^ ~ 1 ^ / sin(q^ cosh( \/l - un)dq j co 8 (qa)da 

sinh ( vT^ q) 




jt^h 


( 22 ) 


o + 1 - u 


is obtained. 

On account of 


cos(qa)da 

+ 1 - li 


• ^ -q\/i^ 

2 \/l - \i ® 


after the introduction of a new variable of integration finally 
becomes 


m„ 




Jth\/l 



sinf 


\ cosh(qT))dq 


Ji w o 


\/l - [ij - 1 


(23) 


The axial additional velocity now reads 

too 

“e / 

u* = 


nh'^d - u) 


q| \ qdq 

C0S{ —==. \cosh(qq) 


.VTt 


1 ^> 


2 q . 1 


( 24 ) 
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The increased velocity at the wall « 0, tj « l) is 

introduced again. The correction velocity at the point ( g = t) = o) 
is then 

,Uo*=ju (25) 

where u is the increased velocity measured at the wall, 

5* DEPENnENCy OF THE DIPOLE MOMENT ON BODI VOLUME AND MACH NUMBER 


The dipole moments and % introduced in equations (3) 
and (19} should not he set in relation to the volume of the body 
that the flow passes any more. Since those potentials only contain 
a single parameter, the volume of the body is the most suitable 
quantity, in fact, for the definition of this parameter. The 
flow past the body could be Introduced as a series of dipoles that 
has set to work in its interior. Each individual dipole signifies 
a certain displacement of the outer flow, which is obtained most 
easily with the aid of the flow function. Therefore, the relation 
between the potential and the flow function must be set up, first 
of all. In its exact form it reads for the case of rotational 
symmetry 

£L M 

Po 

p 3$ (26) 

pQ 


The equations ( 26) are not linear on account of the dependency ' 
between the density p and the velocity, however, they can be 
linearized into the following form:® 


^\}/ 

6 r 



( 26a) 


— : — ^ ^ 

The author is obliged to Dr, Ing, B, Gothert for pointing this out 


The approximate form. i 26 a,)'- ie eqtdval^t to the diffei'ential 
equation for ^ obtained from (l) -fof' 'f the case of 'rotational' 
Bjiumetry and a corresponding eq,^^f fo^ ^ flow function ■■ 

of a dipole in a vmiform flow reads, therefore, in accord with (26a) 




(1 - 


4rt 


\/x2 + (1 - ia)r2 ^ 


(27) 


By Betting equal to 0 the contour of the hody past which 
the stream flows is obtained and frcm this its volume /t. It is 


T = 


2!r 

3U 


(28) 


so that in this case the mcanent is, therefore, independent of the 
Mach number. For the two-dimensional case the relations cox-responding 
to the system (26a) may be written down readily. Fi-om the flow 
fvinction satisfying them 


\|/ = Uy - 


(1 - u)me y ■ 

x^ + (1 - |x)y^ 


(29) 


the volume of the body past which the stream flows (volume within 
surface of the contour past which the stream flows) is obtained as 


J = 


2 U 


(30) 


Fi'om this is obtained the fact that the dipole moment in the 
two-dimensional case is dependent on the Mach number. This result 
is in eiccord with the so-called Prahdtl*s rule (reference 3) - The. 
objection could be raised against this consideration that it 
investigates the flow past a body taking as a basis an individual 
dipole de facto, which does not satisfy the Prandtl condition of 
slenderness. It might, therefore, have been more acceptable to 
represent the body possibly by assuming a distribution of sources 
and sinks along its axis. Now if this is done, then in the extreme 
case of a very slender body admittedly the same dependency of 
the product of source strength by source-eink distance on the 
Mach number is obtained as that for the dipole moments, in, ... _ . . 
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equations (28) and (30), while on the other hand the maierical 
factors change; they hecone equal to 1 in hoth cases, that is 


"" U 

jhq Vl - hi 


(28a) 

(30a) 


The dipole mament for the case of rotational symmetry calculated 
from (28) woul.d then he, accordingly, 50 percent'^ and that for 
the two-dimensional case according to (30) fully 100 percent 
larger than that from the second consideration. Since the bodies 
that appear pi'acticaJ. as models are slender, as a rule, the 
advantage belongs rightly to the second consideration in every case. 

Therefore, introducing (28a) and (30a) into (l5) or (24), now, 
the following is obtained 


u* =■ 


Ut 


t 


CO 


2rt%3 \/i - [i 


/ cos 



l%l(<l) 

Il(<l) 


dq 


(31) 


for the case of rotational symmetry and 



(32) 


for the two-dimensional one. At the position of the body, therefore, 
for y =s p = 0 or g = 11 = 0, the following relations are obtained 

„ * _ 0-1268 

^ It U J _ 0.1309 TJ 3 

^ 2k \/i - u ^ \/ 1 - [I 3 Jj2 


(31a) 

(32a) 
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It is noteworthy that the factor in front of the integral in both 
cases (equations ( 31 ) (3^)) shews tho same dependency on 

the Mach ntaabei’. 


6- BOf.mASE ANGLE CORKECTIOK HJ THE CLOSED TUNNEL AND IN TEE OPEN JET 


In tho two-dimensional case the circulatory flow furnishes 
no conti'ibution to the angle-of-attack correction factor at the 
position of the bocly. On that account only the three-dimensional 
problem in the tunnel of circular cross section is handled in the 
following. For this the action of the model cma bo approximated 
by a horseshoe vortex of infiiiite3.y small span. If instead of 
the velocity potential $ of such a vortex its acceleration 
potential cp is introduced, certain further advantages result, 
in particular, the possibility of keeping the method of solution 
applied up till now. 

The linearized relation between 0 and reads, (reference 4) • 


0 = 


_1 

U 



( 33 ) 


For that very reason (p is also a solution of the differential 
equation (l) . The horseshoe vortex of infinitely small span 
correspond to tho accelez’ation potential of a dipole with its 
axis in tho direction of the 2 -axis (fig. 3); this potential 
reads 


bPU (1 - ^) 2 


cp =■ 


hit \/x^ + (1 - p) ( 




^ Z 




( 31 ^) 


The appropriate velocity potential can bo ascertained frem 
this with ( 33 K It roads 


bi2 


$ 


l(.rt(y‘' 2^) 


i/? + (1 - u) (y'^ + z“) 


+ 1 


( 35 ) 


For given F the Mach number exerts no influence on the 
flow, this holds as much in the plane of the wir.g (x = O) as 
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also Infinitely far ‘behind the wing . At the senw tice it 

is soon that $ and all cross caaponents of the velocity at an 
infinite distance have doiible the value compared to that 
at the position of the wing. The tunnel correction factor at an 
infinite distance la valid, therefore, at the position of the wing, 
too, if it is multiplied "by one-half. . The additional potential of the 
flow ccraing about through the action of the Jet "boundary, at an 
infinite distance, is 



(36) 


in which the upper sign holds for the closed tunnel and the lower 
sign for the open Jet. 

Now the general three -dimensional problem is to be treated. 
At this point an additional potential v* is introduced which 
allows the boundai’y conditions at the edge of the Jet to be 
satisfied. It is easily seen that the boundary conditions(2) 
are also valid 'for the acceleration potential: 


™ (cp + fp*) = 0 
or 


(37) 


Tho boundary condition for the open Jet is obtained as 


q) -i- (p * s 0 


(38) 


The courses of calculation for the open Jot and closed tunnel 
ruin off very much alike. It corresponds, moreover, step by step, 
to the method described in sections 3 and h. An abbi-evieted exposition 
will do here, therefore, in which only the closed tunnel is taken 
up, fii'*st of ell. 

On account of {3h) tho following is applied 


ip* ~ cos 4 P(r) X(x) 



Ik 
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The general solution reads, after making use of (9). 


p = cos ’MCqP + / Cj^ cos 


'll 


kn? 

\ 




(39) 


where la the modified Bessel’s function of the first type and 

first order. 

The coefficients cjf are determined from the hoTHidary condi- 
tions (37)# which are satisfied only for j|| < X. If the 

limit is taken, the final solution is obtained which after 

application of the integral rerresentations for the modified Bessel’s 
functions of the second type takes the form 


-* 

P 


hru cos 1; 
2it^b 2 \/l - p. 



ql %(q) - ql!^(q) 

Ii (qp) 

1 - M Il*(q) 


dq 


(ko) 


Frcm this with the aid of 
is obteinod 


(33) the additional upwash component 


t/* = - 


bP 


I'lr- 


J^OO 


q%(q) - qK (q) 

^ do / cos (qo) dq (4l) 


ii’(q) 


The corresponding upwash component in the open Jet is 


W* = 


bT 


kn^B^- 


'■yi— P d(j cos (qo) 


<0 


I/O 


q\(q) 

I^Cq)" 


dq 


(k 2 ) 


The results (4l) and (42) confirm the observation already made, 
heretofore, that there is no effect due to compressibility at the 
position of the wing and at an infinite distance. In the remainder 
the same additional upwash prevails at a position i behind the 
wing as would be present in incompressible flow at the position 
i 

— ' Since the amount of the correction velocity increases 

\/l - P 

monotonically with increasing distance behind the wing in the open 
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Jet and likewise increases least in the closed tunnel within a 
range comparable to the tunnel radius, the ccmpresslhle flow, 
therefore, has on absolutely larger correction factor. 


7. NUMERICAL RESULTS 


In the following the results of so/eral numerical calculations 
shall be compiled and discussed in detail. The axial velocity 
u* for the case of rotational syaanetry is beet calcinated from 
formula (15) or better still ( 31 ) • For this purpose the integral 


Fi = 



cos(<iD lo(ap) 




dq. 


(*^ 3 ) 


is evaluated numerically by Simpson’s rule ‘^see table l) . 


TABLE 1.' VALUES FOE F. (SEE ( 43 )-) 



p 



0 

0.25 

0.5 

0.75 

1.0 

0 

.25 

.5 

.75 

1.0 

1.5 

0.1268 

.1197 

.1056 

.0853 

.0652 

.0345 

0.1298 

0.1399 

0.1604 

0.1996 

.1877 

.1242 

.0710 

.0409 

.0199 













Next, figure k presents the variation of the additional 
velocity u* along the tunnel radius in the plane x = 0 . Since 
the Mach number in this case only appears in the factor in front 

of the integral, it is sufficient to plot only \/l - |i ^ . 

It is seen that the additional velocity toward the tunnel e“'ge 
takes on possibly 60 percent more. The assumption '■>;» an additional 
velocity V compare reference 5 ), constant over the cross section 
does not prove correct, therefore. 
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u* 

Figure 5 shows the variation of ~ ~ along the tunnel 
axis (r = O) . 

For the velocity at the tunnel wall, frcm a rational point 
of view, not u*, hut the quantity u increased hy the 
displaceinent flow, is plotted for it is certainly this increased 
velocity u which is accessible for direct measurement. The 
var'iation of u as a function of x appears in figure 6. For 
the two-dimensional case (equation (32); it is necessary to 
evaltiate the integral 


1 

Fp =- 
jt 


;os(q|,) cosh(qil) 


The numerical values obtained by Simpson* s “rule are in table 2. 

Figures 7 and 8 show the variation of the additional velocity 
u* along the y-axis (x = O) and along the x-axis (y = O) . 

Figure 9 gives the induced velocity' at the tunnel wall . 

The downwash correction factor for the closed tunnel should 
be represented by means of the upwash w* according to eqxxation (4l) • 
In integrating with respect to a the unsuitable integral can be 
avoided by using the following relation in accord with (36) 


U o 


V . . q%2(-1> *• 

cos(qo) — dq = K 


Then 


(l + Kg) 


where 


Kr, = - 
^ Jt 


\/i-u 


da / cos(qa) 
<J 0 


- q%(q) 
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This fmction has heen tabulated for ■■ ■ -irr rr = 0 ... 5 

Vl - |i 

in which the integral has again been evaluated by Simpson* s rule. 
(See table 3.) 

The curve for the variation of k(j. is shown in figure 10. 
The upwash for the open jet is 


br 

V* ( 1 + %) 

hjiR^ 


by which 




^^(ci) 

IpCl) 


da 


(47) 


(48) 


Table 4 contains the numerical values. 

The ctrves are presented in figure 10. A comparison with 
the variation calculated by I, Lotz (reference l), for u = 0 
and a wing of finitG'wtng span shows good agreement in the case 
of the open jet, on the other hand somewhat larger deviations for 
the closed tunnel, without assigning a reason for this different 
behaviour. On the other hand the variation of both curves of figure 10 
agree very well with the results calculated by Tani and Taima, 
(x-eference 5) using the Bxirgers method. 


SUMMARY 


The problem of the effect of the limitation of the jet on 
the flow past a model is handled by pi'oceeding from the Prandtl 
linearization of the differential equation of the compressible 
medixM,. The disturbance which the model causes near the wall, 
at the same time, is repx’esented, approximately, by a dipole or 
horseshoe vortex. The boundary-value px’oblem ax>ising in this, 
at the limit of the jet is solved exactly to learn tlie additional 
flow due to the effect of the edge of the stream. The solutions 
are evaluated nuraei'*ically, to the extent that they are of interest. 

Translated by Dave Peingold 
National Advisory Committee 
for Aeronautics 
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Figs. 1,2 
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Figure 1.- Designations in the case of rotational symmetry, (x-axis is the 

axis of symmetry). 



Figure 2.- Designations in the two-dimensional case. 
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Figs. 4,5 


Figure 4.- 


Figure 5. 



Axial additional velocity in the plane X = 0* Case of 
rotational symmetry. 





Axial additional velocity along the tunnel axis. Case of 

rotational symmetry. 



Figs. 6,7 
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Lg. 10 



Figure 10.- Downwash correction factors for closed and open tunnels 


